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ABSTRACT 

We present evidence for the existence of an IRAC excess in the spectral 
energy distribution (SED) of 5 galaxies at 0.6 < z < 0.9 and 1 galaxy at z — 1.7. 
These 6 galaxies, located in the Great Observatories Origins Deep Survey field 
(GOODS-N), are star forming since they present strong 6.2, 7.7, and, 11.3 /im 
polycyclic aromatic hydrocarbon (PAH) lines in their Spitzer IRS mid-infrared 
spectra. We use a library of templates computed with PEGASE.2 to fit their 
multiwavelength photometry and derive their stellar continuum. Subtraction of 
the stellar continuum enables us to detect in 5 galaxies a significant excess in 
the IRAC band pass where the 3.3 /im PAH is expected (i.e IRAC 5.8 /im for 
the range of redshifts considered here). We then assess if the physical origin 
of the IRAC excess is due to an obscured active galactic nucleus (AGN) or 
warm dust emission. For one galaxy evidence of an obscured AGN is found, 
while the remaining four do not exhibit any significant AGN activity. Possible 
contamination by warm dust continuum of unknown origin as found in the 
Galactic diffuse emission is discussed. The properties of such a continuum 
would have to be different from the local Universe to explain the measured 
IRAC excess, but we cannot definitively rule out this possibility until its origin 
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is understood. Assuming that the IRAC excess is dominated by the 3.3 fim 
PAH feature, we find good agreement with the observed 11.3 /mi PAH line 
flux arising from the same C-H bending and stretching modes, consistent with 
model expectations. Finally, the IRAC excess appears to be correlated with the 
star-formation rate in the galaxies. Hence it could provide a powerful diagnostic 
for measuring dusty star formation in z > 3 galaxies once the mid-infrared 
spectroscopic capabilities of the James Webb Space Telescope become available. 



Subject headings: Galaxies: evolution - Infrared: galaxies - Galaxies: starburst 



Introduction 



Measuring the star formation history of galaxies as a function of redshift enables the 
build up of stellar mass in the Universe to be constrained. Although there are systematic 
uncertainties between different star formation tracers, results from different studies (UV, 
IR) seem to converge on a flat or a gradual decrease in the star formation rate (SFR) be- 
tween z = 3 and z = 1 followed by a steep er decline between z — 1 and z = (e.g. 
Schiminovich et al.lboOol : Icharv fc ElbAoOlh . 



At high redshift (z > 3), evolution of the SFR still remains uncertain, primarily due to 
the poorly determined dust extinction correction. At these redshifts the SFR in galaxies is es- 



timat ed using the UV luminosity which is strongly affected by dust exti nction (ISteide 



1999) . The extinction corr ection can be calculated using the UV slope (IMeurer et al. 



et al. 



1999 



Adelberger fc Steidell l2000l ) but this technique is affected by several limitations. The UV 
slope can be influenced by the presence of an evolved stellar population and therefore over- 
estimates the extinction correction. Moreover it is well known that local galaxie s harboring 
strong dusty star formation are opaque to UV radiation (e.g. iBuat et al.ll2005l ). This im- 
plies that the UV luminosity may not be a reliable tracer of the SFR in certain galaxies 
and that mid and far-infrared tracers which are correlated with dust emission are required. 
The strongest dust emission features, such as those arising between 6 — 12 /im in polycyclic 
aromatic hydrocarbon (PAH) molecules, are redshifted out of the mid-infrared passband at 
z > 3. As a result, even deep 24 /zm observations are insensitive to dusty star formation 
at these redshifts. Hence, a complete understanding of dust correction will require deep 
observations in the far-infrared/submillimeter regime or require accurate calibration of near- 
infrared tracers of dusty star formation. 

The James Webb Space Telescope (JWST) will obtain spectra between 5 — 27/im and 
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thereby detect the redshifted 3.3/im PAH feature from galaxies at z > 3. The success of previ 



ous s t udies using monochromatic PAH luminos ities as dusty SFR tracers (e.g. lChary &: Elbaz 



200ll ; iTakeuchi et al.ll2005l ; iBrandl et al.ll2006l ) motivates the calibration of this PAH emis- 



sion feature as a SFR indicator. 

Some studies have already investigated the possibility of using the 3.3 jum PAH feature 
as a SFR trac er and have revealed the presence of a correlation betw een L 3 , 3fim and SFR in 
local galaxies (jMouri et al.lll99d ; llmanishill2002ulmanishi et al.ll2006l ). However the existence 
of this 3.3/im PAH feature at higher redshifts needs to be confirmed and its line strength 
calibrated to the true SFR of galaxies. 

In this paper we analyze the multi-wavelength properties of 6 galaxies at z > 0.5 and 
assess the evidence for the presence of the 3.3/im PAH feature. The galaxies display strong 
6.2,7.7, and 11.3 yum PAH emission in their Spitzer mid-infrared spectra and are hence ex- 
pected to also show 3.3/im PAH emission. 

The sa mple of galaxies, situate d in the Great Observatories Origin Deep Survey-North 
(GOODS-N, iDickinson et al.l l2003al ) field, has been observed in the optical by the Hubble 
Space Telescope (HST) and in the infrared by the Spitzer Space Telescope. The multiwave- 
length photometry is fitted with sp ectral energy distributions (SEDs ) from the stellar popu- 
lation synthesis model PEGASE.2 (jFioc &: Rocca-Volmerange!ll997l ). The residual emission 
in the IRAC passband, where the 3.3/im PAH signature should be present, is analyzed to 
determine its origin. We first investigate the possibility that this excess originates from 
an obscured active galactic nucleus (AGN) using the SED of NGC 1068. Then we assess 
the possibility that this excess originates from warm dust emission i.e the 3.3 /im PAH line 
emission and/or a warm dust continuum. To test the hypothesis that the IRAC excess is 
dominated by the 3.3 /im PAH line flux we compare the derived line flux to the 11.3 /im PAH 
flux measured in the IRS spectr a. Since models p redict that the 3.3/im and 11.3/im PAH 
lines originate from C-H modes (ILi fc Draindl200ll ). their line fluxes should be correlated. 

The layout of the paper is as follows: The sample is presented in Section [21 the SED 
fits with PEGASE.2 and the determination of the IRAC excess is presented in Section [3j 
Section H] discusses the origin of the IRAC excess which at these wavelengths can be due to 
hot dust emission from an obscured AGN (Section 14.11) . free-free and/or recombination line 
emission (Section 14. 2p or finally to the 3.3 /im PAH broad emission line and/or a warm dust 
continuum (Section 14.31 and |474"|) . Our conclusions are summarized in Section 



Throughout this paper we will use a cosmology with Hq = 71 km s 1 Mpc 



0.73M 



M 



0.27. 
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2. Sample selection 



The 3.3 fim and the 11.3 \xm PAH lin es arise from stretching and out-of-plane bend- 



ing modes of the C-H bond, respectively (IDuley &: Williamsl Il98ll ; iLeger fc Pugetl Il984j ; 
Li &: Drainei l200ll ). In order to obtain the best constraints on the 3.3 /xm PAH emission 



we consider bright 24/im sources in the GOODS-N field which have spectral coverage of the 
11.3 jum PAH feature. We selected a sample of 21 galaxies with mid-infrared spectra from 
Infrared Spectrograph (IRS; Houck et al. 2004) observations in Spitzer programs GO2-20456 
(PI: Chary) and 262 (PI: Helou). Data reduction includes cleaning rogue pixel s, removing 



laten t charge build-up, removing the sky and averaging the 2D files together (IPope et al. 



20081 ). Spectral extraction was performed using a 2 pixel window in SPICE and data were 
calibrated using the same extraction window on a standard star spectru m. For more detail s 
on the observations and data reduction for the data from GO2-20456, see lPope et al.l (120081 ). 

The data from PID 262 were taken in spectral mapping mode. As for the GO2-20456 
data, the latent charge build up was fit and rogue pixels cleaned. A spectral cube was 



generated using the CUBISM tool (e.g. ISmith et al. 



20071 ). We ensured that we did not 



include the slit loss correction factor (SLCF) while making the cube because the sources in 
the field of view are point sources while the SLCF is used to correct for light diffracted into 
the slit due to a source which is extended wider than the spectrograph slit. The spectral 
cube was generated on a 2.5" spatial grid. The spectrum for each source was extracted 
in a 3x3" square aperture around the source. The flux density in the spectrum was then 
cross-calibrated with the 16 and 24 /zm photometry available for the entire GOODS field. 

We know that the 3.3/im PAH signature, present in the IRAC band, will contribute only 
a few percent of the total IRAC photometry which is dominated by stellar photospheric emis- 
sion. Therefore, the stellar continuum in these 21 objects needs to be accurately estimated. 
We use some of the deepest data currently available in GOODS-N, includi ng 3.6, 4.5, 5.8, 8. 
fim photometry obtained with the Spitzer Infrared Array Camera (IRAC. lFazio et al.ll2004l ). 
and the F435W(S), F606W(U), F775W(z ), F850LP(z) photome try obtained with the HST 
Advanced Camera for Surveys (ACS) (see iDickinson et al.ll2003al ). In order to constrain the 
slope of the stellar continuum in the IRAC passbands which are unaffected by the rising 
dust continuum, we need to apply a redshift cut corresponding to z > 0.5. Moreover the 
redshifting of the 3.3/im PAH feature out of the IRAC passbands implies a second redshift 
cut at z < 1.8. These constraints reduced the initial sample to 11 objects. 

Optical and NIR photometry must be unaffected by neighboring sources since they 
could contaminate the IRAC flux densities by a larger amount than the PAH signature we 
are trying to detect. Using a visual analysis of optical images we rejected every source with 
a neighbor closer than 2 IRAC Full Width Half Max (i.e ~ 3"). This reduces our sub-sample 
of 11 down to only 6 isolated galaxies (see Figured]). Of these 6, only MIPS 3419 and MIPS 
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5581 are from the spectral mapping observations in PID 262. Multiwavelength photometry 
of these galaxies are shown in Table [TJ The optical photometry are MAG_AUTO values 
from the public GOODS-N catalogs. The infrared photometry is measured in 4" diameter 
apertures, with appropriate aperture corrections for the wings of the point spread function. 
Photometric uncertainties are negligibly small (except for MIPS 3419), since t he galaxies are 
very bright, and are do minated by a 5% systematic calibration uncertainty (ISirianni et al. 



20051 : iReach et al J 1200a ). 

To determine the total infrared luminosity of these galaxies, we include th e photometric 



constr aints available from Spitzer 16 / im, 24 /an, and 70 /tm observations from lTeplitz et al. 
( 120061 ). Chary et al. (in prep), and iFrayer et al.l (120061 ) respectively. In order to test the 
presence of an AGN, we also consider Very Large Arra y (VLA) 1.4 GHz obse rvations (Mor- 
rison et al. in prep), and Chandra X-ray observations ( Alexander et al. 2003 1"*! 

All galaxies have measured spectroscopic redshifts. For MIPS 4, MIPS 6, and MIPS 



5581 , redshifts are taken from the Team Keck Trea sury Redshift Surve y (TKRS. IWirth et al. 



2004j), for MIPS 5 and MIPS 7 redshifts are from ICowie et all (120041 ) while for MIPS 3419 



we derived a spectroscopic redshift from its IRS spectra using the 9.7 /im silicate absorption 
feature. Five galaxies (the ones with ground-based spectroscopic redshift, i.e MIPS 4, 5, 6, 7 
and 5581) are at 0.64 < z < 0.84, while the sixth, with the IRS redshift (i.e MIPS 3419), is 
at z — 1.70. The redshift for MIPS 3419 was measured by aligning a local galaxy template 
(Mrk231) whose mid-IR spectral shape is similar to the extracted spectrum and redshifting 
it till it matches the observed spectrum. Since the 9.7/un feature is quite shallow in MIPS 
3419, the range of plausible redshifts in agreement with the data is 1.6 < z < 1.9. This 
is consistent with a visual analysis of the extracted spectrum which shows the flux density 
dropping to zero at ~25/im corresponding to the 9.7/xm silicate feature at z ~ 1.6. The 
galaxies studied in this paper are star forming galaxies since their Spitzer mid-infrared spec- 
tra display strong 7.7 and 11.3 //m PAH emission. 4 out of 6 sources which have adequate 
wavelength coverage in their mid-infrared spectrum even display strong 6.2 /im emission. 
We fit their IRS spectra together with the ir broadband emission at 1 6, 24, 70 /im, and 1.4 
GHz (see Table ED with the SED library of Icharv fc Elbazl (bpOlh (see IPope et aDhoOsl . for 
details). We note that no aperture corrections are necessary since both the spectra and 
broadband observations integrate the entire galaxy. The inferred bolometric luminosities are 
shown in Tableland indicate that they are luminous infrared galaxies (LIRGs), i.e with 
L IR (8 - 1000 /im) > 10 11 L Q and SFR > 17 Moyr" 1 (using the SFR-L IR conversion law). 



1 For a detailed description of all the ancillary data existing in the GOODS fields we refer to the GOODS 
public webpage ( |http: / /www.stsci.edu/ science/goods/| 
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3. Data Analysis 



We estimate the stellar continuum in the IRAC passbands by fitting stellar popula- 
tion synthesis models to the multi-band photometry. The library of stella r emission used as 
templates is computed with PEGASE.2 (IFioc fc Rocca-Volmerangelll997l ). Using the range 
of star formation, infall, and wind histories described in Table __] and available on the Le 
PHARE websitcl (Arnouts et al., in prep), we construct an atlas of templates spanning 
galaxies of all Hubble types. 



For each object, original PEGASE.2 templates are shifted to the spectroscopic red- 
shift of the source and convolved with the transmission of each observed filter. We apply an 
age constraint to ensure that galaxies do not correspond to a template with an age greater 
than the age of the Unive rse at its redshift. F inally we include extinction as a free parameter 
following a Calzetti law (jCalzetti et al.lll994j ). 

The templates are fit to all the BViz photometry and some but not all of the IRAC 
passbands. The fit of the stellar continuum excludes the IRAC bands where the 3.3 /im PAH 
signature may be present and also excludes longer wavelengths which are influenced by the 
rising dust continuum. Hence, for galaxies situated at 0.5 < z < 1, the IRAC 5.8 and 8.0 
/im photometry is excluded from the fit (i.e MIPS 4, MIPS 5, MIPS 6, MIPS 7, MIPS 5581) 
and for the galaxy situated at z = 1.7 (i.e MIPS 3419), only the IRAC 8.0 fim photometry 
is excluded from the fit. 

Applying all these constraints, we calculate the best fit using a \ 2 minimization tech- 
nique. For each object the best fit yields the stellar continuum corresponding to a specific 
combination of 4 parameters: type of galaxy (E, Sa, Sb. . . ), age, normalization (i.e stellar 
mass), and extinction (see Table H]). The result of the individual fitting is presented in Figure 
|2j For all objects the residual difference between the observed photometry and the best fit 
is below 15%. 



Once the best fit is obtained a residual IRAC flux can be derived by computing the 
difference between the IRAC photometry and the stellar continuum (see Tabled]). The re- 
sulting IRAC excess is dependent on the stellar continuum used to fit the data and thus 
relies on the specific value of our 4 parameters (Type of galaxy, Ages, Mass, Extinction). 
To assess this dependence and to estimate the error associated to each IRAC excess that 
we measure, we have undertaken a Monte-Carlo approach. For each source, we randomly 
vary the photometry in each observed band by up to ±3 abend where (Jband is the photometric 
uncertainty in that band. Hence these new values for the photometry are still consistent with 



2 http: //www. oamp.fr/arnouts/LE_PHARE. html 
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the observations. We recompute the best fit SED to these new data values and re-derive the 
inferred IRAC excess. For each object, this procedure is repeated 100 times. The error is 
then defined as the standard deviation of these 100 values. The result of the error estimation 
is presented in Table [5J These error bars account for both the photometric error bars and 
the range of stellar continuum models which fit the observed data. 

We note that the Maraston (2005) models for stellar emission cannot be used for this 
analysis because the empirical spectra of thermally pulsing AGB stars, which are a key fea- 
ture of the models, do not extend beyond rest-frame 2.5 fim. Usage of these models would 
result in an inferred line flux systematically higher by a factor ~2 than the one derived with 
PEGASE.2, due to this break at 2.5 /im. 

The results, shown in Tabled indicate that for all but one galaxy (MIPS 3419), we have 
detected a significant excess in the IRAC 5.8 /mi passband. We note that the absence of an 
IRAC excess for MIPS 3419 could be due to a wrong redshift determination. Indeed as dis- 
cussed in Section [2] the redshift of this source was only derived from its IRS spectrum. In the 
next section, we discuss the physical origin of the IRAC excess for the remaining five galaxies. 



4. Discussion on the origin of the IRAC excess 



The IRAC excess can originate from four different components: (i) hot dust emission 
from an obscured AGN, (ii) free-free and hydrogen recombination line emission from ionized 
gas, ( iii) a pos sible continuum observed in the diffuse m edium of our G alaxy (IFlagey et al. 



20061 ; iLul 12004 ) and in some local star forming galaxies (ILu et al.ll2003l ) and, finally to (iv) 
the 3.3/im PAH feature itself. In the following we discuss the contribution of each of these 
components to the measured IRAC excess. 



4.1. Obscured AGN 



A possible origin for the IRAC excess found in 5 galaxies of our sample could be hot 
dust emission from an obscured AGN. We assess this possibility by studying the X-ray, radio, 
and optical properties of these galaxies. 

Using Chandra observations in GOODS-N Alexander et al.ll2003h we find that three 
galaxies of our sample (MIPS 4, 6, and 7) are detected in the soft X-rays (0.5-2 KeV) and 
only one in the hard band (2-8 KeV) (MIPS 6; see Table [6]). The photon index T = 1.06 
of MIPS 6 indicates the presence of an o bscured AGN. Moreo ver, the fact that the resolved 
fraction of the X-ray background (CXB, IWorsley et al.ll2004l ) decreases with increasing en- 



- 8- 



ergy suggests the existence of a population of o bscured AGN which might be missed in even 
the deepest X-ray surveys (IBarger et al.l 120071 ). To address the presence of such AGNs we 
compare the SED of the galaxies in our sample with the SED of the Compton thick AGN 
in NGC 1068. The choice of NGC 1068 has been made since its SED follows a F\ oc \~ a 
power-law with a = 2.25 in the NIR. This is a typical spectral index for obscured AGN 
rtRisaliti et aJbood ). 

The SED for the nucleus of NGC 1068 was deriv e d from the high spatial resolution 
ground-based photometry compiled b y iGalliano et al.l (120031 ). the ISOPHOT-S spectrum 
presented by iRigopoulou et al.l (119991 ) . and the IRAS broadband photometry. The IRAS 
large beam photometry at 12, 25, 60, and 100 /on w as fit with a far- i nfrare d curve compris- 
ing of multiple temperature dust components as in IChary &: Elbazl (120011 ). This synthetic 
spectrum was nor malized at 12 nm to the I SOPHOT-S spectrum of the nucleus which extends 
between 6 — 12 /im (IRigopoulou et al.lll999l ). The ISOPHOT-S spectrum is clearly dominated 
by continuum emission from the AGN. However, the flux density appears to systematically 
exceed the high resolution ground-based photometry by a factor of 2. It is unlikely that this 
is due to extended emission in the vicinity of the nucleus which enters the ISOPHOT-S aper- 
ture since the spectrum is clearly dominated by hot dust emission. It is more likely due to 
flux calibration uncertainties associated with the ISOPHOT-S spectrum. After dividing the 
integrated ISOPHOT+IRAS spectrum by a factor of 2, bringing all the data in agreement, 
we extend the spectrum, using the ground-based photometry, down to 1 /im. 

Assuming that the excess observed in the IRAC passband is only due to hot dust 
from an obscured AGN, we calculate the normalization factor for the NGC 1068 SED corre- 
sponding to this excess. Then subtracting the AGN contribution from the multi-wavelength 
photometry and using the fitting procedure described in section [3j we fit the revised photom- 
etry with the PEGASE.2 population synthesis model. The result of these fits are presented 
in Figure [31 

For 3 objects, MIPS 4, 5, and 5581, these fits yield a predicted flux in the IRAC 8 /im 
passband which exceeds the observed photometry by a factor 6.2, 4.7, and 18.6 a respec- 
tively. These over-estimations of the IRAC 8 /im photometry prove that MIPS 4, 5 and 5581 
cannot harbor such an obscured AGN since their IRAC color (IRAC5.8/IRAC8.0) do not 
follow the typical colors of hot dust emission in NGC 1068. Moreover to match the IRAC 
8 /xm photometry the obscured AGN, if present in these galaxies, would have to be at least 
a factor of 1.3, 1.25, and 2.0 less luminous respectively. Then these obscured AGN would 
not be luminous enough to explain the 5.8 /im IRAC excess. On the contrary, for MIPS 6 
and 7, the renormalized NGC 1068 SED results in an IRAC 8 /tm flux which is consistent 
with observations and thus suggests that hot dust from an obscured AGN could explain the 
IRAC excess. One should note that these conclusions are dependent on the NIR spectral 
index of NGC 1068 (i.e a = 2.25). To assess this dependence we also attempted to fit the 
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IRAC excess using the SED of Mrk 231 (jArmus et al.l 120071 ) which has a low spectral index 
of a = 0.25. We find that the AGN fit exceeds the observed 8 /mi IRAC flux for MIPS 4, 5, 
and 5581 by a value of 4.1, 3.0, and 9.9 a respectively. Although these values are a factor of 
two smaller than the ones derived for NGC 1068 they also weaken the possibility that hot 
dust emission can explain the IRAC excess. 



We now evaluate if the observed X-ray emission from these sources is consistent with 
this simple model. We assume that the X-ray emission is the sum of starburst and AGN 
contributions, and compare with the Chandra observations (see Table EJ). Since our galaxies 
show strong PAH features, the bolometric luminosity calculated in Section 2, is dominated 
by star-formation. We calculate the st arburst X-ray contribution using the relationship 
between SFR and L IR (IKennicuttl Il998l ) and the relationship between SFR and L 5 _ 8keV 
feauer et aDl2002f ): 

SFRlMvyr- 1 } = 4.5 x lO^Lr^ergs- 1 ] (1) 
SFR^yr' 1 ] = 1.7 x lO^L^vfergs- 1 ] (2) 



Using T = 1.9 as the typical starburst photon index we estimate the soft and hard 
X-ray rest-frame luminosity (£go/t/hord) P re dicted for the starburst. The observed flux in 



each band ( f 



Bauer et al 



observed \ 
ard) 



is then derived by applying a i^-correction give by Equation (2) of 



r observed r - 
J soft/hard L 6 



^cm" 2 ! 



T restframe r _n 
L soft/ Hardin S ] 



4vrd 2 



1+z 



i r-2 



(3) 



where di is the luminosity dist ance. For the AG N contribution, we adopt L 2 _i k e v = 
2.8 x lO^ergcm^s" 1 and V = 1.2 Jorie et alJl2003h for NGC 1068. From this, we predict 
the emission of our objects in the soft and hard X-ray bands (cf Figure [3]) using the NGC 
1068 normalization factor calculated from the IRAC excess and the i^-correction given by 
Eq. m 



For all galaxies the total (AGN + Starburst contribution) soft X-ray emission, dom- 
inated by the starburst contribution (L^*y£ st / ~ 10), is in agreement, within the 
error, with the observed soft X-ray flux by a factor of 1.3, 1.7, 1.8, and 1.0 for MIPS 4, 5, 6, 
and 7 respectively and under the Chandra threshold for MIPS 5581 (see Table [6]). 

For MIPS 4, 5, 7, and 5581, the total hard X-ray prediction is under the Chandra 
threshold and thus in agreement with observations. For MIPS 6, our hard X-ray prediction 
(1.48 x 10~ 16 erg s _1 cm -2 ) is a factor of 1.7 below the observation (2.54 x 10~ 16 erg s _1 cm -2 ) 
but still, within the error, in agreement. This suggests that the obscured AGN harbored by 
this galaxy is significantly more luminous than for NGC 1068 in the X-rays. 

In summary, the X-ray analysis reveals the presence of an obscured AGN in MIPS 6 
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but is not conclusive for MIPS 4, 5, 7, and 5581 since either the presence or absence of an 
obscured AGN would be in agreement with the data. The X-ray predictions using the Mrk 
231 SED are similar to those derived using the SED of NGC 1068. Hence, our conclusions 
are independent of the AGN SED adopted to fit the data. 



Using the Chary & Elbaz library which follows the local radio/far- infrared correla- 
tion we predict L 1AG n z using the 24 /im luminosity (see Table [2]). For MIPS 4, 5, and, 5581 
we find L^ r ^^ d / Ll b l GUz ~ 0.8, and for MIPS 6 and MIPS 7, which may be powered by an 
AGN, we find L\ r lo^ / L^i GHz ~ 1.2 and 0.8 respectively. This agreement between predic- 
tions and observations suggests that the infrared emission of all these galaxies is dominated 
by star-formation. Finally, assuming a spectral index (a in F u oc v~ a ) of 0.8 we calculate 
for each galaxy ^i^ghT™ 6 ( see Table ED- These values are then compared to 5 x 10 23 W.Hz _1 
which is the typical ra dio luminosity u s ed to distinguish b etween a starburst and an AGN 
population at z ~ 0.7 (jYun et al.ll200ll ; ICowie et al.l I2004T ) . For each galaxy in our sample, 
^x^GHz" 16 ^ s below 5 x 10 23 W.Hz~ 1 which is consistent with the assumption that our galaxies 
are star-formation dominated systems. 



Finally, we analyze the optical spectra of three galaxies (MIPS 4, 6, and 5581) available 
in the TKRS database (jWirth et al.l I2004J ) . No clear evidence of high ionization lines like 
[Ne III] or [Ne V] or large [OIII]/iJg ratio are found arguing against the possibility that the 
IRAC excess of these three galaxies is dominated from AGN activity. However, we note that 
this analysis includes MIPS 6 which shows clear evidence of an AGN in its X-ray emission. 



The conclusion of the SED fitting, X-ray, optical and radio diagnostics on the presence 
of an obscured AGN in these galaxies can be summarized as follows: 

- For MIPS 4, 5, and, 5581 the disagreement between the IRAC 8 /im expected for a 
prototypical obscured AGN and the observations suggests that these galaxies do not harbor 
an AGN. Moreover for MIPS 4 and 5581, this conclusion is also supported by the agreement 
found between their radio and mid-infrared luminosities as well as the absence of AGN sig- 
nature in their optical spectra. 

- For MIPS 7, the IRAC excess could be reproduced by a buried AGN but we found 
no other evidence for an AGN neither from its optical nor its X-ray properties. 

- Finally, MIPS 6 presents marginal evidence for the presence of an obscured AGN in 
its X-ray emission (1.7 times larger than typically expected for star- formation). Its radio 
luminosity and optical spectrum are consistent with a star forming dominated system. In 
the following, we will discuss the case of this galaxy separately since we cannot rule out its 
contamination by an obscured AGN. 
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4.2. Free-free and gas lines 



The total contribution of free-free and recombination line emission have been estimated 



bv iFlagev et all (120061 ) to be about 1% and 3% in the IRAC 3.6 /iin (i.e the residual IRAC 
5.8 /im emission for MIPS 4, 5, 6, 7 and 5581 or the residual IRAC 8.0 /im for MIPS 3419) 
and in the IRAC 4.8 /im (i.e the residual IRAC 8.0 /im for MIPS 4, 5, 6, 7 and, 5581) 
channels respectively. Even in the extreme case where these contributions would reach 3% 
and 11% (IFlagev et al.ll2006l ). the free- free and line emission contribution to the IRAC excess 
would be negligible. Furthermore if the IRAC excess was due to free-free emission, it would 
be accompanied by a comparably high value at shorter wavelengths (i.e in the 4.5 /im IRAC 
channel) due to the blue spectrum of free-free emission. The flux density in the 3.6 /im and 
4.5 /im IRAC channels in MIPS 4 ,5, 6, 7 and 5581 are consistent with being dominated by 
stellar emission (see Figure [2]) excluding this possibility. 



4.3. Warm dust continuum 



A contin uum underlying the 3.3/im P AH feature has been detected in the Ga lactic dif- 



fuse m edium (IFlagey et al.ll2006l ; lLull2004j ) and in some local star forming galaxies (ILu et al. 



20031 ). There, the NIR continuum, well fitted by a modified black body with temperature 
spanning the range 700 — 1500 K and a X~ 2 emissivity law, contributes to about 70% of the 
IRAC 3.6 /im flux. They also conclude that the IRAC 4.5 /im channel is totally dominated 
by this NIR continuum. 

Assuming that the IRAC excess, that we measure here in distant galaxies, is due 
to such a continuum, the fP redicted ( 5. 8 fim) / fP redicted (8.0 /im) color for MIPS 4, 5, 6, 7 
and, 5581 would span the range [0.97 - 2.07]. We can assume these color predictions 
to be lower limits since they assume no contribution of the 3.3/im PAH feature in the 
IRAC 5.8 /im passband. The observed ratios after subtracting the stellar contribution 
(i.e [/° b ™ d (5.8/im)-^ PBG ^ are 0.64, 

0.64, 0.45, 0.32 and 1.08 for MIPS 4, 5, 6, 7 and 5581 respectively. Although, one galaxy 
is marginally consistent with this warm dust continuum (MIPS 5581), the remaining four 
correspond to a modified black body tempe rature of 550 K, hen ce colder than the temper- 
ature range of the NIR continuum found in IFlagev et al.l (120061 ). As a result it is unlikely, 
unless the shape of the continuum were different from anything known, that the IRAC 5.8 
fim passband is dominated by the same NIR continuum in MIPS 4, 5, 6 and, 7. 

The physical origin of such a high temperature continuum remains unclear and we 
cannot totally rule out that it reaches a different temperature range in distant highly star 
forming galaxies. Near-infrared spectroscopy (L-band, 34 /im) of a sample of 24 local ac- 
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tively star forming galaxies (jlmanishi & Dudleyll2000l ) does exhibit some contribution from 
a hot continuum which could, in principle, serves as a better reference for the present sam- 
ple of distant star forming galaxies. We find that t heir continuum can be reproduced by 
a colder modified black-body than the one found in iFlagey et al.l (120061 ) since it reaches a 
temperature of T~ 400 K. This colder continuum, which is more consistent with the IRAC 
temperature of the present sample, would contribute for 60% of the IRAC excess. However, 
this sample consists only of the nuclei of nearby ultraluminous infrared galaxies (ULIRGs) 
obtained using the Subaru IRCS near-infrared spectrograph. These galaxies were chosen as 
candidate AGNs, the goal of the authors being to look for buried AGNs, which they claim 
to find with various intensities in all objects. Hence the warm dust continuum observed in 
their sample can be highly contaminated by AGN activity. 



4.4. 3.3 fim PAH feature 



Assuming that the IRAC excess is due to the 3.3 jum PAH feature, we show that the 
inferred line fluxes are consistent with expectations from dust models. To infer the 3.3 fim 
PAH line flux from the IRAC excess, we adopt a Gaussian profile for the PAH line. Then 
we convolve this profile with the IRAC filter curve to obtain the normalization factor for the 
Gaussian profile corresponding to the IRAC excess. Finally the true flux of the 3.3/im PAH 
line is obtained by integrating over the normalized Gaussian profile. The estimated value 
of the 3.3 fim PAH line flux calculated for each galaxy is shown in Table [5j We note that 
the line flux inferred by adopting a Drude profile is only 1.12 times the value derived for a 
Gaussian profile. Therefore, our conclusion does not depend on this choice. We also note 
that the extinction derived from the stellar population fits to the photometry (see Table 
|4|) reaches an average value of v4(3.3 fim) fa 0.14. Hence, the extinction correction to the 
inferred 3.3 fim PAH lines fluxes is smaller than the uncertainty in the flux and are therefore 
negligible. 

Models predict that the 3.3 fim and th e 11.3 fim lines originate in the same C-H bond 
( )Li fc Drainei l200ll ; iDuley fc Williams! 119811 ) which results in a linear correlation between 
these two lines. Hence, finding a correlation between the inferred 3.3 fim line flux and 
11.3 fim line flux, measured from the IRS spectra, would support the assumption of a negli- 
gible continuum. 



Using the iDraine fc Lil (120071 ) templates and a Gaussian profile for the PAH features, 
we find Ls^^m = aLn^^ m with a = 0.4 ± 0.2 in the models. However these templates do 
no t represent t h e extr eme case where the PAH are totally ionized or neutral. Using Table 1 
of iLi fc Drainei (120011 ). we find a = 0.3 or 1.3 for the ionized and neutral PAH respectively. 
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The 11.3/xm line flux of each galaxy has been measured in its IRS spectrum using 
the ISAP (jffigdon et al.l l2004h component of SMART (see Table EJ. The line fit is done 
assuming a Gaussian profile for the PAH line and a constant continuum across the line. The 
uncertainty in the line flux is dominated by the signal-to-noise and wavelength coverage in 
the IRS spectrum which makes the derivation of the continuum slope difficult. In the case 
of MIPS 7, where the wavelength coverage does not encompass the entire 11.3/ xm PAH line. 



we si mulate the line flux uncertainty by using spectra of low redshift galaxies (lArmus et al. 



20071 ) and applying a similar wavelength cutoff. We find that the maximum uncertainty in 



the derivation of the line flux is due to the unknown shape of the underlying continuum and 
is at most 40% for MIPS7. Figure H] shows the comparison between the 3.3 /im and 11.3 /im 
PAH line flux and prediction from the models. 

The 3.3 /im line flux value for MIPS 3419 is under 2 a hence, it has been considered as 
an upper limit. For comparison, we plot observations of two star-formation dominated local 
galaxies (Arp 220 and Mrk 273) galaxies and 2 AGN (Mrk 231, IRAS 05189-2425) taken from 



Imanishi fc Dudley! (120001 ) and lArmus et al.l (120071 ) for L 3 . 3Atm and Lu ^ m respectively. The 



classification of Arp 220 and Mrk 273 as star-formation dominated galaxies an d of Mrk 231 



and IR AS 05189-2425 as AGN are based on the mid-infrared diagram sh own inlArmus et al. 
(2007J, see their Figure 8). We note th at the 3.3 /im line fl ux derived by llmanishi fc Dudley 



( l2000l ) includes only the nucleus while lArmus et al.l (120071 ) measure the 11.3 /im feature in a 
larger aperture which encompasses most of the galaxy. We apply an aperture correction to 
the 3.3 /xm line flux which is the ratio of the ~ 3 — 4 /im continuum of the whole galaxy ob- 
tained using broad band observations reported in th e NASA/IPAC Extr a galac tic Database 
(NED) and the 3 — 4 /im continuum in the spectra of llmanishi fc Dudley! (|2000| ). This aper- 
ture correction relies on the assumption that the nucleus and the whole galaxy have the 
same spectrum which may not be the case for AGN-dominated nuclei. In fact, we note that 
this aperture corrections would provide only a lower limit to L 3 3fim due to the low 3.3 /im 
equivalent widths observed in obscured AGN. 

The values of L 3 3fJim inferred for 5 out of 6 ga l axies (MIPS 4, 5, 6, 7 and, 3419) are con- 
sistent with the one predicted by the iLi &: Draind (120011 ) model from the measured Ln.3 Mm . 
This agreement supports the hypothesis that the IRAC excess observed in these galaxies is 
due to the 3.3 /im PAH feature. This true also for MIPS 6 and MIPS 7 suggesting that their 
potential AGN contribution (see Section 14.11) does not dominate the IRAC flux excess. For 
MIPS 5581, the inferred 3.3 /im PAH line is over-estimated by a factor of ~ 1.5 — 4 with 
respect to model expectations from the Ln, 3fJirn . This can be explained by the presence of 
a significant NIR continuum as also revealed by the color diagnostic presented in Section 
14.31 The assumption that the NIR continuum is negligible compared to the line leads to an 
overestimation of the inferred 3.3 /im PAH line. 

The presence of strong PAH features in the IRS spectra of each galaxy of the sample, ex- 
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cept for MIPS 3419, indicates that they must be star-formation d ominated. The L3 ^ m/LfR 
ratio in local star- forming galaxies has been found to be ~ 10~ 3 (jMouri et al.lll990l ) with a 
factor of 2 - 3 scatter. For MIPS 4, 5, 6, 7, and 5581 we find L^ m /L IR ~ 3 x 10" 3 which 
is, within the uncertainties of the relation, consistent with the local ratio. For MIPS 3419, 
the upper limit found for the 3.3 fim PAH line flux is still consistent with the local ratio of 

L3.3fj, m / Lir. 

The agreement of the L3.3 Mm /In.3 Mm ratio with standard dust models and of the 
£3.3^111/ L>m ratio with local observations favors the hypothesis that the IRAC excess is ef- 
fectively due to the 3.3 fim PAH emission. We also note that even in the case where the 
IRAC excess would be due to the combined e mission o f a con tinuum and the 3.3 /im PAH, 
the study of nearby star-forming galaxies by iLu et al.l (120031 ) showed that both emissions 
probably originate from the same carriers. Indeed, they find that the 7.7 //m PAH equiva- 
lent width d o es no t vary in their sample of galaxies. This trend was not corroborated by 
Flagey et al.l (120061) but t he lat ter study remains limited to the diffuse ISM of the Milky 
Way. In light of lLu et~aL ( 2003 ) results the information carried by the 3.3 fim PAH line flux 
alone or by the continuum and the 3.3 /im PAH complex would be equivalent. 



5. Conclusion 

We have investigated the presence of an excess in the observed IRAC flux densities 
through analysis of the multi-wavelength photometry and spectroscopy of 6 galaxies in the 
redshift range 0.5 < z < 1.8. Using a fit to the optical and NIR photometry, we determined 
the stellar continuum of each source with PEGASE.2. The difference between the stellar 
continuum and the IRAC flux in the passband expected to harbor the 3.3 /im PAH feature 
is computed and analyzed. For 5 galaxies we found a significant IRAC excess. For the other 
galaxy, we were only able to determine an upper limit. 

We investigated the possibility that the measured IRAC excess could be explained by 
the presence of an obscured AGN using the SED of NGC 1068. For 1 galaxy (MIPS 6) we 
find evidence for a possible contamination by an obscured AGN. For 4 galaxies (MIPS 4, 5, 
7, and 5581), no evidence was found that the IRAC excess could result from the presence of 
an obscured AGN. 

We have investigated the origin of the IRAC excess as potentially due to the presence 



of th e 3.3 jLtm PAH fe ature and/or a warm dust continuum. Local observations ( jFlagey et al. 



20061 ; ILu et al.l 120031 ) find evidence for a warm dust continuum of unknown origin which 
could dominate the IRAC excess. The NIR colors of the IRAC excess cannot be reproduced 
by such a continuum unless one assumes a lower temperature. However, the 3.3 fim PAH 
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line flux inferred directly from the IRAC excess is consistent with standard dust models. 
Indeed, we find that the inferred 3.3 /mi PAH line flux is compatible with the 11.3 /mi PAH 
line flux measured in the IRS spectrum. This suggests that the IRAC excess is effectively 
due to the 3.3 /tm PAH feature. 

However, we note that even if the IRAC excess would arise from the combination of the 
3.3 (Urn PA H line flux and a warm dust continuum such as that found in local star-forming 
galaxies by lLu et al.l (120031 ). both originate from the same carriers, as indicated by a constant 
PAH/continuum flux ratio. In that well as in the case where the 3.3 /zm PAH line 

flux alone dominates the IRAC excess, we conclude that this emission provides a powerful 
diagnostic for measuring dusty star formation rates in galaxies at z > 3 using the JWST. 
We find that uL^f l^^/LjR is nearby constant for all 5 galaxies, i.e (1.0 ±0.3) x 10~ 3 with 
Ljpt spanning the range [1.5 x 10 11 — 8 x 10 11 ]. 
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Fig. 1. — From left to right and from top to bottom each panel shows optical and NIR images of MIPS4, 
MIPS5, MIPS6, MIPS7, MIPS3419, and MIPS5581 respectively. For each galaxy, we show on the left a color 
composite BVz image from HST and on the right IRAC 5.8 /jm from Spitzer. Each image spans a 10" x 10" 
area. 
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Fig. 2. — Determination of the stellar continuum for 6 galaxies in our sample. Optical and near-infrared 
photometry used to fit the stellar continuum are shown by red squares, near-infrared photometry excluded 
from the fit are shown by red stars, and the IRS spectra by a continuous line (all IRS spectra have been 
smoothed to the same spectral resolution of R ~ 45 at 15 fim, using a Gaussian kernel). The stellar continuum 
computed by PEGASE.2 is shown with the red dotted line. The convolution of the continuum with each 
passband is shown as the open circle (only wavelengths blueward of the shaded region have been used to 
determine the fits to the stellar continuum). The dashed vertical line represents the observed wavelength 
at which the redshifted 3.3 fim PAH features would be present and the shaded area represents the bandpass 
of the IRAC filter in which this PAH features falls. Thin vertical lines show the location of the 6.2, 7.7, 
and 11.3 /zm PAH lines. The bottom of each panel shows the ratio of the residual error, defined as the 
difference between the observations and the final fit, and the photometric uncertainty in each passband. 
Errors obtained for the passbands used to fit the stellar continuum are shown as black squares while black 
stars represent errors obtained for the passbands excluded from the fit. The difference at red wavelengths is 
most likely due to the 3.3 /im PAH and the increasing contribution from the VSG continuum. 
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Fig. 3. — Determination of the stellar continuum for 5 galaxies from which we have removed an AGN 
contribution proportional to the IRAC excess calculated in section [3] Lines and symbols are the same as 
in Figure [2l The stellar continuum computed by PEGASE.2 and the AGN contribution are shown as dot 
and dash lines respectively (as before, only wavelengths blueward of the shaded region have been used to 
determine the stellar continuum). The final fit, sum of the stellar and the AGN contributions, convolved 
through each band pass is shown as open circle. The X-ray contribution of the normalized AGN is shown in 
each panel in ergs" 1 cm" 2 . For sources 4, 5 and 5581, the IRAC excess is inconsistent with being due to an 
obscured AGN (see text for details). 
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Fig. 4. — Comparison between the inferred 3.3 /im PAH line flux and the 11.3/im PAH line flux measured 



from the IRS spectra. Uncertainties for L 



3.3 fim 



are computed using a Monte Carlo approach (see text 



for details). The arrow corresponds to a 2a limit. The das hed lines represent L3 
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with 



a = 0.3 or 1.3 for the ionized and neutral PAH respectively (|Li fc Draind 120011) . For comparison, we plot 
the PAH luminosities of 2 loca l starburst galaxies (A rp 220 and Mrk 27 3) and 2 local AGN (Mrk 231, IRAS 
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Table 1. Optical and near- infrared photometry of our sample 



Objects 



Id 


z 


RA 


Dec 


B 


V 


i 


z 


3.6 fim 


4.5 £tm 


5.8 (an 


8.0 






(J2000) 


(/«Jy) 


(MJy) 


(**Jy) 


(**Jy) 


(MJy) 


(^Jy) 


(i«Jy) 


(^Jy) 


MIPS4 


0.638 


189.01355 


62.18634 


3.9 ±0.2 


11.9 ±0.5 


27.2 ± 1.2 


35.8 ± 1.6 


129.9 ±6.5 


89.1 ±4.4 


107.9 ±5.4 


98.7 ±4.9 


MIPS5 


0.641 


189.39383 


62.28978 


1.0 ±0.06 


5.0 ±0.2 


15.4 ±0.7 


24.3 ± 1.0 


167.9 ±8.4 


113.9 ±5.7 


124.9 ±6.2 


100.9 ±5.0 


MIPS6 


0.639 


189.09367 


62.26231 


5.3 ±0.2 


12.1 ±0.5 


24.5 ± 1.1 


30.4 ± 1.3 


67.3 ±3.3 


46.5 ±2.3 


56.4 ±2.8 


65.8 ± 3.3 


MIPS7 


0.792 


189.23306 


62.13559 


0.3 ±0.06 


2.0 ±0.09 


8.5 ±0.3 


13.5 ±0.6 


112.9 ± 5.6 


80.8 ±4.0 


77.9 ±3.9 


80.9 ±4.0 


MIPS3419 


1.70 


189.17568 


62.28963 


0.02 ±0.01 


0.06 ±0.01 


0.35 ±0.03 


0.55 ±0.04 


33.5 ± 1.6 


38.5 ± 1.9 


34.3 ± 1.7 


23.5 ± 1.2 


MIPS5581 


0.839 


189.28491 


62.25418 


2.4 ±0.1 


3.9 ±0.2 


7.9 ±0.3 


9.9 ±0.4 


27.3 ± 1.3 


18.7 ±0.9 


20.0 ± 1.0 


14.7 ± 0.8 



Note. — Catalogs have been cross-correlated using a matching radius of 0.5 
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Table 2. Infrared and radio flux densities of the sample 



Id object 


(fJy) 


524 jim 

(fJy) 


570 /im 
(mjy) 


qobs 

"1.4GHz 

(fJy) 


^predicted 
J 1.4GHz 

(MJy) 


T rest 
-■"1.4GHz 
(xl0 23 WHz _1 ) 


Log(L IR ) 

(L ) 


MIPS 4 


777 ± 20 


1221 ± 12 


11.0 ±0.66 


161 ± 12 


126 


2.4 


11.90 ±0.05 


MIPS 5 


575 ± 11 


750 ±7 


5.6 ±0.67 


91 ± 14 


79 


1.4 


11.75 ± 0.04 


MIPS 6 


398 ±5 


721 ± 7 


11.0 ± 0.66 


64 ± 16 


75 


1.0 


11.70 ±0.04 


MIPS 7 


582 ± 10 


832 ± 8 


14.0 ±0.7 


104 ± 11 


81 


2.7 


11.89 ± 0.05 


MIPS 3419 


54 ± 7 


113 ± 6 


< 3.0 


< 25 


19 


< 4.0 


11.76 ± 0.10 


MIPS 5581 


194 ±5 


201 ± 6 


< 3.0 


16 ±5 


17 


0.6 


11.22 ± 0.09 



Table 3. PEGASE.2 template parameters 



Type v Infall (t c ) Gal Winds 

(Myr) (Myr) (Myr) 



E 100 100 3000 

SO 500 100 5000 

Sa 1500 500 

Sb 2500 1000 

Sbc 5000 1000 

Sc 10000 2000 

Sd 20000 2000 

Irr 20000 5000 



Note. — PEGASE.2 scenarios used as 
template parameters. SFR = u —1 X 
Mgas and gas infall is simulated as /(t) = 
exp( - t/t c ) ^ ^ j s effectively the ratio be- 
tween the star formation time scale and the 
star formation efficiency. The Initial Mass 
Function used in ou r scenarios is taken from 
iRana fc Basul l|l992h . 
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Table 4. Fit parameters 



Id object 


Type 


Age 


Stellar Mass 


E(B-V) 


9 

x 2 






(Gyr) 


(M ) 






MIPS 4 


Sbc 


7 


1.8 x 10 11 


0.18 


5.15 


MIPS 5 


Sa 


6 


2.8 x 10 11 


0.72 


9.75 


MIPS 6 


Sbc 


5 


7.3 x 10 10 


0.03 


8.68 


MIPS 7 


E 


3 


2.3 x 10 11 


0.88 


15.75 


MIPS 3419 


Sd 


3 


2.6 x 10 11 


2.73 


6.48 


MIPS 5581 


Sb 


3 


3.4 x 10 10 


0.03 


22.26 


Note. — x 2 


— 2^i=i 


(*i-Si) 2 
of 


- where N is 


the number 


of pass- 


bands fit to the 


model. 











Table 5. Inferred PAH properties 



Id object IRAC excess 


3.3/im line flux 


11.3/um line flux 


EW r e st frame 


SFR 


L3.3 nm/LiR 


(xlO _2S crg s _1 cm _2 Hz _1 ) 


(xlO _22 Wcm -2 ) 


(xlO~ 22 Wcm _2 ) 


H 


(M Yr" 1 ) 


(X10~ 3 ) 


MIPS 4 3.7 ±0.7 


4.92 ± 1.45 


4.67 ±0.58 


227 


137.9 


2.7 ±0.8 


MIPS 5 3.2 ±0.8 


3.66 ±0.75 


5.87 ±0.68 


58 


97.8 


2.81 ±0.8 


MIPS 6 2.0 ±0.3 


2.75 ±0.30 


3.69 ±0.69 


119 


86.3 


2.48 ±0.5 


MIPS 7 1.4 ±0.5 


1.34 ±0.62 


4.16 ± 1.66 


45 


136.0 


1.29 ±0.6 


MIPS 3419 0.0 ±0.1 


0.00 ±0.11 


< 0.18 


< 11 


99.5 


< 1.9 


MIPS 5581 0.6 ±0.1 


0.67 ±0.15 


0.17 ±0.9 


106 


28.5 


3.5 ± 1.0 


Note. — 3.3 fj,m line fluxes have been calculated with the original PEGASE.2 fit and errors have been calculated using the 


Monte Carlo approach. The SFR is derived from Lju using Equation [T] 








Table 6. X-Ray properties of the sample 






Id object 


Soft X-ray 


Hard X-ray 


Photon index 








(0.5-2 keV) 


(2-8 keV) 


(T) 








Obs (Model) 


Obs (Model) 










[xl0 _16 crg 


; s — 1 cm -2 ] 








4 


1.42 (1.86) 


< 3.21 (2.33) 


> 1.45 






5 


< 0.76 (1.31) 


< 7.49 (1.68) 








6 


0.65 (1.22) 


2.54 (1.48) 


1.06 






7 


1.06 (1.07) 


< 2.73 (1.32) 


> 1.35 






5581 


< 0.29 (0.24) 


< 1.55 (0.33) 









